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The Nambu-Jona-Lasinio (NJL) - Jet model provides a sound framework for calculating the frag- 
mentation functions in an effective chiral quark theory, where the momentum and isospin sum rules 
are satisfied without the introduction of ad hoc parameters. Earlier studies of the pion fragmenta- 
tion functions using the NJL model within this framework showed qualitative agreement with the 
empirical parameterizations. Here we extend the NJL-Jet model by including the strange quark. 
The corrections to the pion fragmentation functions and corresponding kaon fragmentation func- 
tions are calculated using the elementary quark to quark-meson fragmentation functions from NJL. 
The results for the kaon fragmentation functionsa exhibit a qualitative agreement with the empirical 
parameterizations, while the unfavored strange quark fragmentation to pions is shown to be of the 
same order of magnitude as the unfavored light quark's. The results of these studies are expected 
to provide important guidance for the analysis of a large variety of semi-inclusive data. 
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I. QUARK FRAGMENTATION AND NJL-JET MODEL 

Quark fragmentation functions have long been of interest for analyzing hard scattering reactions [IHZ] • New exper- 
imental efforts to extract distribution and fragmentation functions from deep-inelastic lepton-nucleon, proton-proton 
scattering and e + e~ annihilation data [5HTT] have generated a renewed interest in the subject [T2JIT3]. The analysis 
of the transversity quark distribution functions [5J Q3] and a variety of other semi-inclusive processes [T51 [TFJ also 
critically depend on the knowledge of the fragmentation functions, both unpolarized and polarized [T71 118j . 

The NJL-Jet model of Ref. [19] provides a self-consistent framework for calculations of both quark distributions 
and fragmentation functions in an effective chiral quark theory. The advantage of the model is that there are no ad 
hoc parameters introduced to describe fragmentation functions. The isospin and momentum sum rules are satisfied 
naturally in the product ansatz. In this work we extend the NJL-Jet model by introducing the strange quark, thus 
allowing fragmentation to both pions and kaons. The inclusion of the new channel is shown to bring significant 
corrections to the previously calculated fragmentation functions to pions leading to a better agreement with empirical 
parametrizations. 

II. INCLUDING THE STRANGE QUARK 
A. Strange Quark Constituent Mass and Coupling 

Throughout this article we employ the SU(3) NJL model with a straightforward generalization of the SU(2) 
Lagrangian with no additional free parameters (see for e.g. Refs. [20H22] for detailed reviews of the NJL model) 
and we consider only the pseudoscalar mesons. Thus introducing the strange quark in the NJL-Jet model involves 
calculating the strange quark distribution and fragmentation functions in the NJL framework, requiring the knowledge 
of the quark-meson (strange-light to kaon) coupling constant and the strange quark constituent mass. We calculate 
the quark-meson coupling using the same approach as Ref. [23] • The strange quark mass is chosen to best fit the 
model calculated kaon mass to the experimental value. 

The quark-meson coupling constant is determined from the residue at the pole in the quark-antiquark t-matrix at 
the mass of the meson under consideration. This involves the derivative of the familiar quark-bubble graph: 
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Here Tr denotes the Dirac-trace and the subscripts on the quark propagators denote quarks of different flavor - also 
indicated by q and Q, where the meson of type m under consideration has the flavor structure m = qQ. We evaluate 
Eq. ([I]) in light-cone (LC) coordinates 1 , noting: 
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The II(p, ar) is evaluated using the complex residue theorem and we obtain for p± = 0: 
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Thus for the quark-meson coupling we obtain: 
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The integrals in the above expressions are divergent and require regularization. Here we use the Lepage-Brodsky 
(LB) "invariant mass" cut-off regularization (see Refs. [HI |23j for a detailed description as applied to the NJL-Jet 
model), where A 12 is the maximum invariant mass of the two particles in the loop: 



Mu < A 12 = yjAl + M? + sJk% + M%, (8) 
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where A3 denotes the 3-momentum cutoff, which is fixed in the usual way by reproducing the experimental pion decay 
constant. We use the value obtained in the Ref. [53] of A3 = 0.67 GeV. 

Thus from the requirement of k\ > we obtain the region of x where the integrand is non-zero. 

Using the light constituent quark mass M = 0.3 GeV from the Ref. [23] and the experimental value of kaon mass 
771 k — 0.495 GeV yields a strange constituent quark mass M s = 0.537 GeV and the corresponding quark-kaon coupling 
constant of gxqQ = 3.39. 

B. Quark Distribution and Fragmentation Functions 

The quark distribution function fg(x) has an interpretation as the probability to find a quark of type q with 
momentum fraction x in the hadron (in our case meson) h. The corresponding cut diagram is shown in Fig.JI^,, which 
can be equivalently represented by the Feynman diagram depicted in Fig. [T]d: 



1 We use the following LC convention for Lorentz 4- vectors (a + ,a , aj_), = "^j( a ° ± a 3 ) and ax = (a 1 , a 2 ). 
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FIG. 1: Quark distribution functions. 
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FIG. 2: The model solutions for distribution functions of quarks in mesons at a) model scale Q% — 0.2 GeV 2 and b) evolved at 
next to leading order to scale Q 2 = 4 GeV 2 . 
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where k_ = xp_ and C™ is the corresponding flavor factor given in the Table |TJ Here we used the identity 
S(k)j + S(k) = —dS(k)dk + to relate the distribution function to the function Tl(p, x) and the expression in the 
Eq. (|6| to evaluate the result. 

In our calculations we use the same LB regularization scheme (|8| yielding the non-zero region for f^_(x) as described 

in the previous section. The corresponding plots for , f*f and fjf are depicted in Fig. ^ at model scale 
Qq = 0.2 GeV 2 and evolved at next to leading order (NLO) to scale Q 2 = 4 GeV 2 using the software from Ref. |24) . 
Here the model scale is set by reproducing the behavior, after NLO evolution, of the u quark valence distribution 
function in the pion given in Refs. |251 126) . 

It is easy to check using Eqs. ([3 11) that both isospin and momentum sum rules for the distribution functions are 
satisfied: 



dxf?(x) = -C™, (12) 
1 dxx (f?(x) + f%(xj) = \C™, m = qQ, (13) 

where in the last line we used the relations /™(x) = — x) (which is easy to verify from Eq. |Tl| ) and C™ = 
(see the Table [I]). 
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FIG. 3: Quark fragmentation functions. 

The elementary fragmentation function depicted in Fig. [3] in the frame where the fragmenting quark has kj^ = 
(but non-zero transverse momentum component fcy = —p_\_/z with respect to the direction of the produced hadron) 
can be written as: 

W = -^ mqQ \ J j- )4 Tr[S 1 (kh+S 1 (k)^-p + M 2 ) l5 ] 

x6(k_ -p_/z)2ir5{{p-kf - M|) = -^-J?(x = 1/z) (14) 
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The integration can be trivially done analytically if one assumes a sharp cut-off in the transverse momentum, 
denoted by P]_ ( in the LB regularization scheme P]_ is z-dependent): 

{z) = °l 9 l^z(4lif^ + lo g (l + Pl/B)), (16) 



2 8tt 2 \B/Pl 
where 



A = {(z- 1)M X + M 2 ) 2 , (17) 
B = z(z - 1)M 2 + zM\ + (1 - z)m 2 m . (18) 



In LB regularization the P\ is given by: 



K =z{l-z) U/A£+n& + y/AZ + Mij -(l-z)m z m -zMi. (19) 
The corresponding sum rule for the elementary fragmentation function is as follows: 

f dz d™(z) = 1-Z™. 
Jo 



(20) 



Here Z™ is the residue of the quark propagator of flavor q in the presence of the cloud of the meson type m. It is 
expressed in terms of the renormalized quark self-energy induced by the meson type m loop £™(fc) of Fig. |4] as 

£"(*) = -. C™gl, Q J i£f,7sS«(«: - P )7 5 A,„(p), (22) 
where M q is the mass of the constituent quark of flavor q and A m (p) is the Feynman propagator of the meson type 

TO. 



5 



k ft w k 
> • » #— > 

FIG. 4: (Color Online) Quark Self-Energy. 
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FIG. 5: Quark cascade. 

III. GENERALIZED NJL-JET 

The NJL-Jet model of Ref. [TH] uses a multiplicative ansatz for the total fragmentation function to derive an 
integral equation for the quark cascade for the process depicted in Fig. [5j The derived integral equation for the total 
fragmentation function is: 



D™(z) = <(z) + J2 f -^(-) D%(y), df(z) = d™(l - z)\ m=qQ . (23) 

Q Jz v y 



Here d™(z) — d™(z) / Ylm' — ^™ ) ( nere ^he sum is over all possible meson tpyes m' that the quark of flavor q can 
emit in the elementary splitting process). Then Y) m J d™(z)dz — 1, thus allowing an interpretation as the probability 



of an elementary process. In Eq. (23) the sum is over the flavor of the emitted quark Q and the splitting function 
of quark q into Q with a light-cone momentum fraction z is naturally the same as the splitting function of q into a 
meson m of a flavor composition qQ with a light-cone momentum fraction 1 — z. 



The result in Eq.(23) resembles the integral equation ansatz of Field and Feynman's quark-jet model [U |2"7] . 



Rewriting the above expression helps to elucidate the probabilistic interpretation of the model: 

D™(z)dz = d™(z)dz + J2 d%j)dyD%{-)-. (24) 
Q Jz V V 

Here the left hand side term has the meaning of the probability to create a meson m carrying the light-cone 
momentum fraction z to z + dz of initial quark q. The first term on the right hand side corresponds to the probability 
of creating the meson with light-cone momentum fraction z to z + dz in the first step of the cascade, plus the second 
term corresponding to the creation of the meson further down the quark cascade after a splitting to a quark Q with 
light-cone momentum fraction y. Here the probability of creating the meson m with the light-cone momentum fraction 
z of the initial quark is the probability of creating the same meson with the light-cone momentum fraction z/y in the 
cascade of the quark Q, which is clearly only the case in the Bjorken limit. Thus the model can be generalized in a 



straightforward way by including the strange quark directly in Eq. ( 23 1 . 

We solve the coupled set of integral equations using the elementary fragmentation functions of Eq. (23) for u, 
d and s quark fragmentation to a given meson. The corresponding fragmentation functions for the anti-quarks are 
obtained using charge symmetry starting from the fragmentation functions of quarks to the corresponding anti-meson. 
The comparisons with the phenomenological parametrizations of Ref. [TU] are performed by DGLAP evolution of the 
calculated fragmentation functions from the low-energy model scale of Ql = 0.2 GeV 2 to 4 GeV 2 , at NLO, using the 
software from Ref. [24] , 
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FIG. 6: The total fractions of momenta carried by mesons of type m from the jet of a) u and b) s quarks calculated using the 
numerical solutions for fragmentation functions. 



It is easy to see, using the properties of d™ along with the normalization condition of , that the solutions of the 
integral equations (23) should satisfy both momentum and isospin sum rules. We use a notation: 



and define: 



(f(z)) / f(z)dz, 
Jo 
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(25) 



Using the equations Eq. (23) it follows: 



N = n + n'N. 



(26) 



(27) 



To prove the momentum sum rule iV = 1, that is the fragmenting quark transfers all of its momentum to the 
emitted mesons, it is sufficient to show that n + n' = 1: 



n + n 



' = E ({<{*)) + (^(i - *))) = E = 1 



(28) 



The total fractions of momenta carried by mesons of type m from the jet of u and s quarks, ^zZ?™(z)^) ; q = 
u, s calculated using the numerical solutions for fragmentation functions are shown in Fig. [6j which show that the 
momentum sum rule is satisfied within the numerical precision of the calculation of less than a percent. 

Similarly we define t as the z-component of the isospin and 



At q = Y,t m (D™(z)) , at q = ^2t m (d™(z)) , a't q = (d$(z) 

m m Q 



Using the equations Eq. ( 23 ) it follows 



A = a + a' A. 



(29) 



(30) 



To prove the isospin sum rule A = 1, that is the z-component of the isospin is conserved in the quark fragmentation 
cascade, it is sufficient to show that a + a' = 1: 
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at q + a%= J2 ({th%{z)) + (t Q d?(l-z)))=t q ^(Z?{z))=t q , (31) 

TO 

to = qQ 

where we kept the flavor of the fragmenting quark q constant and the sum is over all possible mesons to that can be 
emitted (for each to the corresponding fragmented quark is denoted Q). 

Our numerical solutions obey these rules within numerical errors of less than a percent, with, for example, 86% of 
the u quark's isospin transferred to pions and 14% to kaons. 



IV. RESULTS AND CONCLUSIONS 



The results for the fragmentation functions of u, d and s quarks to 7r + and K + at the model scale are shown 
in Fig. [7] Though the fragmentation functions of the unfavored strange and light quarks are of the same order of 
magnitude, we can see a notable difference between them, even at the low scale of the model. This is potentially 
significant observation as it has hitherto been necessary in phenomenological analyses to assume that Dl = D* a , 
etc. (e.g. see Ref. [28] ) . 

Further we present our results for the fragmentation functions D™ + and = D? + in Fig. ^ Here the DGLAP 
evolved curve is also compared to the empirical parametrizations of the experimental data of Ref. [10] , evolved to 
the same scale. We can see that the inclusion of strangeness softens the high z region of Z?J + compared to previous 
calculations [T5], thus bringing the curves closer to the phenomenological parametrizations (here we shall note that 
the model fragmentation functions were evolved at NLO in the current article versus the LO evolved solutions of the 
Ref. [TH] )• This is expected as the elementary fragmentation to a kaon is not negligible. In fact the fragmentation to 
K + is about half as likely as that to n + as can be seen from the plots in Fig. [9] The plots show a reasonably good 
agreement with the parametrizations within the large uncertainties of the latter. 

It is clear that for a more complete description of the quark fragmentation both vector meson and nucleon anti- 
nucleon channels need to be included in the calculations. This can be accomplished within the current framework. 
The high z region of fragmentation functions are dominated by "few-step" transitions where the availability of the 
additional fragmentation channels might have a noticeable effect. 

Another limitation of the model is the assumption of the momentum scaling of the probability of hadron creation 
in each step of the decay chain, which is clearly only the case in the Bjorkcn limit. A quark with a finite momentum 
loses energy with each production of a hadron and finally recombines with the remnants of the antiquark jet to form 
the final hadron. A more accurate description of the process requires Monte-Carlo (MC) simulations of the quark 
fragmentation, similar to the studies in Refs. (23[29]j and others. MC simulations would also allow one to access the 
transverse momentum distribution of the produced hadrons, thus becoming relevant for the analysis of a large variety 
of semi-inclusive data. 
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FIG. 8: Pion fragmentation functions for a) favored zDJ (2) and b) unfavored zD^ (z) process . Here the dashed line with 
the error band represents the empirical parametrization of the Ref. [10], the doted line is the solution for the NJL-Jet model 
at the low energy scale Q 2 , = 0.2 GeV 2 and the solid line represents the same solution evol ved at next to leading order to the 
scale Q 2 = 4 GeV 2 . The dash-dotted line is the driving term z d™(z) in integral equation (23 1. 
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FIG. 9: Kaon fragmentation functions for a) favored zD„ (z) and b) unfavored zDff (z) process . Here the dashed line with 
the error band represents the empirical parametrization of the Ref. [TO], the doted line is the solution for the NJL-Jet model 
at the low energy scale Q 2 = 0.2 GeV 2 and the solid line represents the same solution evol ved at next to leading order to the 



scale Q = 4 GeV . The dash-dotted line is the driving term z d™(z) in integral equation (23 1 
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Appendix A: Flavor Factors for Distribution and Splitting Functions 



The flavor factors given in the Table [T] are calculated from the corresponding SU (3) flavor matrices (for details see 
for e.g. Ref. [21]): 
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